Scalar mesons moving in a finite volume and the role of partial wave mixing 
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Phase shifts and resonance parameters can be obtained from finite-volume lattice spectra for 
interacting pairs of particles, moving with nonzero total momentum. We present a simple derivation 
of the method that is subsequently applied to obtain the tt-k and nK phase shifts in the sectors with 
total isospin I — and / = 1/2, respectively. Considering different total momenta, one obtains extra 
data points for a given volume that allow for a very efficient extraction of the resonance parameters 
in the infinite- volume limit. Corrections due to the mixing of partial waves are provided. We expect 
that our results will help to optimize the strategies in lattice simulations, which aim at an accurate 
determination of the scattering and resonance properties. 
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I. INTRODUCTION 

One of the present issues in QCD lattice calculations is 
the determination of the excited hadron spectrum. Many 
efforts arc being devoted to this problem lately jlPr27l]. 
In the volume-dependent spectrum, the "avoided level 
crossing" is usually taken as a signal of a resonance, but 
this criterion has been shown insufficient for resonances 
with a large width [28j - [32| . For resonances with a single 
decay channel, one often uses Liischer's approach to ex- 
tract phase shifts from the discrete energy levels in the 
box [H, [HI- The method has been recently extended 
to multi-channel scattering [3l|, |35T - |37| and to the case 
with three-particle intermediate states [38l . [39j . More- 
over, in Ref. [3^|, a method based on coupled-channel 
Unitary Chiral Perturbation Theory (UChPT), which op- 
erates with the full relativistic two-body propagator, has 
been proposed. In the infinite-volume limit, this method 
is equivalent to Liischer's approach, up to contributions 
(kept in Ref. [32| ) , which are exponentially suppressed in 
this limit. The new method, combining conceptual and 
technical simplicity, can provide a guideline for future 
lattice calculations. 

This method has been extended in Ref. for the 
use in connection with dynamical coupled-channel ap- 
proaches such as the Jiilich model |4l| - [43| . In Ref. [44[ 
the framework of Ref. [12] has been applied for the inter- 
action of the DK and rjD s system, where the D s «o(2317) 
resonance is dynamically generated from the interaction 
of these particles. The case of the k in the Ktt, Krj chan- 
nels is addressed along the lines of Ref. [32| in Ref. [45[ , 
together with the <r(600), K*(892), and p resonances. 
The formalism has also been extended to the case of the 
interaction of unstable particles in Ref. [46| , to the study 
of the DN interaction [47| , and the 7T7t interaction in the 
p channel [48| . Pioneering work for the coupled-channel 
KN, 7rS system and the A(1405) in the finite volume 
has been carried out in Ref. [33] ; the lattice levels for the 



A(1405) quantum numbers were evaluated in Ref. (40[ 
using the Jiilich model, and strategies to determine the 
two A(1405) states from lattice results (c.f. also Ref. 0) 
were discussed in Ref. [i^. Methods to extract matrix 
elements of unstable particles from the finite volume have 
been recently developed in Ref. J5C 

32j is done 



The derivation of Refs. [33j, [34| or Ref. 
for a pair of particles with total zero-momentum. The 
genera lization to a moving frame has been done in [39l . 

In this study we present an easy derivation of the ap- 
proach for moving frames, along the lines of Ref. [33 ] . 
using fully relativistic propagators and arbitrary masses, 
and we apply the method to study the coupled-channel 
scattering of 7T7t and irK in the region of the /o (600) and 
k(800) resonances, respectively. 

The strategy followed here is to use the chiral unitary 
approach to generate synthetic lattice data, which are 
later on analyzed to extract phase shifts and resonance 
properties. We show that the data produced with moving 
frames are very useful to get the infinite-volume proper- 
ties using two box volumes and different total momenta. 
The study done here permits to find optimal strategies, 
concerning which lattice data to use, in order to obtain 
phase shifts and resonance properties with maximum pre- 
cision. 



II. FORMALISM 
A. Particles in a moving frame 

In the chiral unitary approach the scattering matrix in 
coupled channels is given by the Bethe-Salpeter equation 
(BSE) in its factorized form 



T = [1 - VG^V = [V- 1 - G}'\ 



(1) 
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where V = V^' is the matrix for the transition poten- 
tials between the channels and G is a diagonal matrix 
with the i th element, G^' , given by the loop function of 
two propagators, which for two mesons is defined as 
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(2tt) 4 (P - q) 2 -mj+ie q 2 - m\ 



(2) 

where nii are the masses of the two mesons and P the 
total four-momentum of the meson-meson system. The 
factorized form of the BSE implies an on-shcll factoriza- 
tion of the potential V, see Sec. Ill Cl for a discussion. 

V and T in Eq. ([T]) stand for the potential and scatter- 
ing matrix in the momentum space. The normalization 
is such that, for the case of one channel, we have 
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which relates T to the phase shift S and the S'-matrix, 
where E is the total energy in the CM and p the mo- 
mentum, p = \ 1/2 (E 2 ,ml,ml)/(2E). For the two- 
channel case, corresponding relations can be found in 

Refs. m, m. 

The loop function in Eq. ^ needs to be regularized 
and this can be accomplished cither with dimensional 
regularization or with a three-momentum cutoff. The 
equivalence of both methods was shown in Refs. [6lLl62T|. 
In the regularization with a three-momentum cutoff, one 
first performs the q° integration analytically [63|. As a 
result, one gets 



W\«ln 



where 
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We would like to stress that other renormalization 
schemes such as dimensional regularization are, in gen- 
eral, preferable over the cut-off renormalization. How- 
ever, note that for the extraction of the infinite- volume 
limit from lattice levels the cut-off dependence cancels 
as discussed in Sec. Ill Gl The treatment of finite volume 
using dimensional regularization is done in Ref . |44| . 

To obtain the energy levels in the finite box, instead 
of integrating over the momenta of the continuum - with 
q being a continuous variable as in Eq. ^ - one must 
sum over the discrete momenta allowed in a finite box 
of side length L with periodic boundary conditions. We 



then have to replace G by G — diag (G^\ & 2 ^), where 

|«|«!W 



2n 



q = —n, n € 

lJ 



(0) 



This is the procedure followed in [32j. Here and in the 
following, we indicate quantities in the finite volume with 
a tilde, e.g. G ->• G, or T -> f. 

The eigenenergies of the box correspond to energies 
that produce poles in the T scattering matrix in the finite 
volume, 



f =[l-VG]- 1 V =[V~ 1 -G] 



(7) 



i.e. for energies where det(l — VG) = 0. 

In the former discussion the integrals and sums are 
performed in the rest frame of the two interacting parti- 
cles. Yet, since in the infinite volume the G function is 
Lorentz invariant, see Eq. it suffices to evaluate it in 
the two-particle rest frame. In another frame it will take 
the same value as required by Lorentz invariancc. Let 
us call q* = \q*\ the absolute value of the relative three- 
momentum in the rest frame of the two particles, called 
center-of-mass frame (CM) in the following. The absolute 
value of the relative three-momentum in a frame where 
the two-particle system has total momentum (P°,P) is 
called q. The CM energy of the two-particle system will 
be y/s, such that 



= E 2 
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However, for the system moving in the finite volume, 
Lorentz invariance is broken and hence we cannot use 
G evaluated in the CM, because the discretization condi- 
tion in the momenta q{ and q 2 = P — q[ of the particles 
of Eq. ^ must be transformed to the moving frame. We 
must write the boost transformation from q to q* . By ap- 
plying the Lorentz transformation from a moving frame 
with four- momentum P to a frame where the two particle 
system is at rest we find 



<?1.2 = 91,2 



H-l) ^- 2 ' P _ ^1 

J IPI 2 " \Ts 



P. (9) 



Here and in the following, a star indicates a quantity 
defined in the two-particle rest frame. 

Demanding that q\ + g| =0 enforces q± + q 2 = P°. 
We also have the transformation of the energies 



o _ fa*^ + gi.2 ■ p) 

and the condition q® + q 2 = P° imposes 



qf + qf 



(10) 
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which on shell gives 
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This, via Eqs. (J9j) and (flQ|) . provides then the boost for 
the off-shell momenta in the loop, where q is arbitrary 
but the energy is the on-shcll one. Only this prescription 



<72 for two particles of equal mass in the 



ensures q{° 

two-particle rest frame. Since we need the Jacobian of 
this transformation, it is useful to rewrite Eq. ^ in terms 
of the CM energy of the particles and we find 



Ql.2 = <7i,2 + 
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This equation is the one used in [641 ]. Furthermore 
we must substitute / d 3 q* /(2tt) 3 by / d 3 q/(2ir) 3 y/s/P , 
where the factor y/s /P° is the Jacobian of the transfor- 
mation, and then replace the integral by the discrete sum. 
In summary, we must perform the substitution 
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where for both the sum and the integral the limit is \q* \ < 

ymax- 

Note that in order to have both q{ and q^ = P — q[ ful- 
filling the periodic boundary conditions, the momentum 
P must fulfill them, too, and thus we have 
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(16) 



A clarification on the role of the on-shcll reduction used 
in Eq. ([T]) is appropriate. In general, one does not know 
how good this approximation is. However, in the infinite- 
volume limit, the on-shell approximation of the potential 
V derived from a Lagrangian delivers in many cases a 
successful description of the phenomenology, see, e.g., the 
hadronic model we use in this work [62| . However, fully 
covariant off-shell unitarized approaches have been also 
developed [65T4671 ] . Off-shell unitarized approaches in a 
three-dimensional reduction are realized, e.g., in dynam- 
ical coupled-channel models [4ll - |43j . 

While there are a conceptual differences between on- 
and off-shell approaches in the infinite volume limit, it 
should be clearly stated that, whatever approach is cho- 
sen in the infinite volume limit, the levels in the finite 
volume are determined by the on-shcll amplitude up to 
exponentially suppressed effects ~ e~ L M,r . 



infinite volume can be obtained for the energies which 
are eigenvalues of the box by (E = t/s) 

T{E) = (V-\E) - G(E))- 1 = (G(P) G(E))- 1 . 

(17) 

since G{P) = V 1 (E) is the condition for the T matrix 
to have a pole for the finite box. 

Hence we find, in the one-channel case and assuming 
that only S'-wave scattering is present, 
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This derivation is very simple and the results can be 
seen to agree with previous ones [Hll [HI, [55|, [5(| when one 
approximates I(q*) by 
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2E p 2 - {q*Y 



(18) 



where p = X 1 / 2 (E 2 , mf, m\)/{2E) and X(x,y,z) stands 
for the Kallen triangle function. It is seen that, with 
this replacement, our expressions agree with those in the 
Liischer framework. Furthermore, summing the differ- 
ence of Eq. and Eq. (|18[) over the momenta, it is im- 
mediately seen that the finite- volume corrections to this 
quantity are exponentially suppressed. Consequently, the 
present approach is equivalent to Liischer's approach up 
to the exponentially suppressed terms for large volumes. 
Note however that, for moderately large volumes, these 
exponentially suppressed terms can be important numer- 
ically, see the discussion in Ref. [32j . One should realize, 
however, that for values of LM^ where these terms play 
a role, there are many other exponentially suppressed 
corrections. 

We would like to stress that the formal dependence 
on the potential V cancels in Eq. (|TT)) - in other words, 
Eq. (jTTJ) contains T and does not contain V or G indi- 
vidually. This was expected from the beginning, because 
the potential V is not an observable and depends on the 
cut-off chosen for G. For the case of moving frames and 
partial wave mixing, discussed in the following, we will 
find exactly the same behavior: the measured lattice lev- 
els depend only on the T-matrix in the infinite volume 
and and functions G — G, which are cutoff-independent. 



C. Partial wave decomposition in a finite volume 



B. One-channel analysis 

The one-channel problem can be easily solved and is 
very simple, as shown in [32| . The T matrix for the 



To determine the mixing of partial waves, consider first 
the case without boost, i.e. P = 0. At the end of this 
section, the formalism is generalized to moving frames 
and multiple channels. 
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We use the spherical harmonics Y^ m with the normal- 
ization 

/•27T l"K 

/ dcf> / sm6 d0Y em (6, cf>)Y £ 1 m ,(6, cf>) = 5 u >5 mm , , 
Jo Jo 

(19) 

and further define 

yimip) = p e Y em (0, 0) = p e Y £m (p) , 

p = p (sin 9 cos (/>, sin 9 sin (p, cos 9) = pp . (20) 

The potential V in the Bcthe-Salpeter equation (fTJ) is the 
same in a finite and in the infinite volume, V = V . Its 
partial-wave expansion takes the standard form 

V(p, p') = All Y, if) Vi (p, p') y* £m (f ') • (21) 

However, the partial-wave expansion of the scattering 
amplitude T is different in a finite volume, because here 
the rotational symmetry is broken down to cubic sym- 
metry. As a result, 

f(p,p') = 4tt]T yi m (p)t im , e > m '(p,p')yz m '(p')(M) 

£,m 
£'m' 

In the above expressions, v £ and t(_ m ^i m i depend only 
on p 2 = p 2 and (p 1 ) 2 = (p 1 ) 2 . Note that the thresh- 
old behavior of the amplitude is hidden in the function 

ye m (p) ~ P l - 

In the infinite- volume limit, the rotational symmetry is 
restored and the Wigncr-Eckart theorem guarantees that 
t is diagonal both in I and m, 

tlm.l'm' = ti^U'Smm,' ■ (23) 

In contrast, in the finite volume one obtains 

Um,£'m'(P,p') = Ve(p,p') Su'Smm' + 

ft l"m" 

(q,p') (24) 

with I(q) from Eq. ([5]). This relation is obtained straight- 
forwardly by inserting V from Eq. (|21[) and T from 
Eq. 1221) in Eq. (0). 

Next, note that, since vt{p,p') and te m ,i> m '(p,p') de- 
pend only on p 2 ,(p') 2 , the use of the regular summation 
theorem [33| is justified. Further, the argument q in these 
functions [see Eq. (|2"4"|)] can be replaced by the on-shell 
value q on , which is determined from the zero of the en- 
ergy denominator in Eq. © . Finally, using Eq. (|2(J|) , one 
gets 

Ttm.i'm'iPiP') = Vi(p,p')8 U 'Sm m ' 

+ ]T Ve(p,q o nGirn,e>> m ^qnT£» m »,e' m >(q on ,p'), 

t"m" 

(25) 



where 

Te m ,i'm'(p,p') = p e te m ,£' m '{p,p'){p'Y , 

V,{p,p')=p l v l {p,p'){p') i , (26) 

and 

G tm , em , (q° a ) = ^ E (^) + F /™(9) *W (q) ■ 

ft 

(27) 

The factor (q/q° n ) e+t can be replaced by 1, when £ + £' 
is even, and by q/q on otherwise [38|], i.e. 

(q/q on Y +e ' -+ (q/q on ) k , 

k = 0, 1 for £ + £' = even, odd . (28) 

Note that only even £ + £' lead to non-zero contributions 
for P = 0, as well as in case of equal-mass particle scat- 
tering. 

Note also that if the factor (q/q on ) is neglected one 
obtains expressions that are different from the origi- 
nal Liischer approach [33|, [H[ not by exponentially sup- 
pressed terms, but terms suppressed as 1/L 4 (38l |. In any 
case, we have checked numerically that effects for the re- 
sults of this study, coming from the (q/q on ) factor, are 
very small for the considered realistic box sizes. 

An important remark is in order. At first glance, it 
seems that there is an ambiguity in the choice of k in 
Eq. I|28p. Note, however, that this problem arose be- 
cause one insisted on the on-shell prescription in the in- 
finite volume limit. The above choice ensures that our 
finite-volume expressions are compatible with Luscher's 
approach [33|, Hi[ up to exponentially suppressed terms, 
and in the infinite-volume limit they are also compatible 
with the on-shell prescription. In particular, the quantity 
G in the infinite-volume limit is replaced by 

f d 3 q 

G£m,e>m'(q° n ) > Sti'Smm' / , n N o I (q) 

L-roo J {2-K)- i 

= o~ee>S mm > G , (29) 

with G from Eq. @ (channel index omitted here). The 
Bethe-Salpeter equation in the infinite- volume limit takes 
the simple form given in Eq. (TT]), i.e. 

T e (p,p>) = V e (p,p') + V t (p, q° n ) G(q oa ) T e (q° n ,p>) (30) 

where the arguments are quoted explicitly. Finally, recall 
that the expressions in Eqs. ([2T| and (|2"9")) arc defined with 
an implicit momentum cutoff at q = <7 max . 

D. Partial wave mixing with boost and multiple 
channels 

In Sec. Ill CI we have considered the partial-wave ex- 
pansion for a zero total momentum, P = 0, of the two 
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particles. Nothing changes conceptually if we consider 
moving frames instead, and the formalism of Sec. Ill Al 
can be applied. In this case, the quantity G is given by 



4-7T 



E 



E 



P \q on * 



xY; m (r)Y e , m „(q*)I(q*) 



(31) 



Here, q* is given by Eq. (fT5|). E = yfs is from Eq. ©, 
fc from Eq. (J2HJ) , and Ge m ,e"m" is obviously the general- 
ization of G from Eq. (jT5J) to higher partial waves. Note 
that, in case of particles with different masses, the states 
with even values of £ and odd values of I' can mix in a 
moving frame. In this case, as discussed above, we choose 
k = 1. 

In the large- L limit, up to the terms exponentially sup- 
pressed in L, the quantity Ge m ,£"m" can be expressed in 
terms of Liischer's zeta-functions. Namely, 



G 



lm,l"m" 



8u"5mm"G — — 



871-yi 



M 



ira.t" m 1 



(32) 



where Mt m ,t"m" is given, e.g., by Eq. (39) of Ref. [73j. 
It is a linear combination of the Luscher zeta-functions 
in the moving frame. 

The discrete levels in a finite volume emerge at the en- 
ergies where the determinant of the linear equation (|25|) 
vanishes, 

det (Su'Smm' - W\ 0&m,<'m'(0) = ■ ( 33 ) 

Here, G is a matrix with a row (column) index given by 
all combinations of £, m (£', m'). 

Finally, all above formulae refer to the single-channel 
case. In case of multiple coupled channels, both the po- 
tential V and the quantity G should be considered as 
matrices in channel space. In that case, V, which is di- 
agonal in the space of partial waves, obtains additional 

(ii) 

channel indices Vg — > V t . The quantity G, that has 
non-diagonal elements in the space of partial waves due to 
the mixing, becomes a diagonal matrix in channel space, 

Gi m ,e'rn> -> Gy^ e , m ,, c.f. also Eq. (J2J). 

In the following, along with UChPT, we shall use 
the Inverse Amplitude Method (IAM) in the version 
of Ref. [62| in the finite volume, as has been done al- 
ready in Ref. (H for the case P = 0. The IAM 



method exploits unitarity, which in one channel states 
that Im T- 1 = p/{8nE) (see Eq. ©). Denoting by 
and the second- and fourth-order chiral po- 
tentials, and noting that the amplitude at lowest order, 
ym = T^ 2 \ has no imaginary part, a dispersion rela- 
tion is made for the function (V^) 2 /T, where the imag- 
inary part is thus known analytically. It leads to a sim- 
ple relationship T IAM = (tW) 2 /(TW -T^), where T®, 
and TW are the amplitudes at lowest order and next- 
to- lowest order, respectively j6^|. We can then recast 
these results in terms of Eq. |T]), redefining a potential 
y[4] = T [4] _ y[2\ GV [2\_ The equivalent potential to be 

used in Eq. |T]) or (f3"3"| to obtain the results of the IAM 
amplitude is then given by 



V, 



IAM 



i-v^ty^)- 1 ) v [2] 



(34) 



Here, = (V^)^ and V® = (yM)fa') are matrices 
in channel space. For V^ 4 \ only the polynomial terms 
are considered like in Ref. [62j. The method derived 
here could be also extended to the loop calculations of 

Refs. Ei-tzi. 



For illustration, an explicit example for the emerging 
structure is quoted. With V ( M = (Vj AM )W> , we con- 
sider the mixing of S- and P-waves in the coupled chan- 
nels wK and r/K. Then, V and G in Eq. (|33p are given 
by 
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and 



6 



/ r (1) 
"00,00 





Mi) 

"00,1-1 





Mi) 

"00,10 





Mi) 

"00,11 


\ 





Mi) 

"00, 00 





M*) 

"00,1-1 





/=V(2) 
"00,10 





M?) 

"00,11 


M 1 ) 

"1-1, 00 





M 1 ) 

"1-1,1-1 





Mi) 

"1-1,10 





Mi) 
"1-1,11 








M i 

"1-1,00 





M A) 

"1-1,1-1 





M A) 

"1-1,10 





M i 
"i-i,ii 


Mi) 

"10,00 





Mi) 

"10,1-1 





Mi) 

"10,10 





Mi) 

"10,11 








M2) 

"10,00 





M*) 

"10,1-1 





M2) 
"10,10 





M2) 
"10,11 


Mi) 

"11,00 





Mi) 
"11,1-1 





Mi) 

"11,10 





Mi) 
"ii,ii 





V o 


"11,00 





"ii,i-i 





M2) 
"11,10 





M2) , 
^n, n / 



E. Symmetries of Eq. ([33]) 

As mentioned above, rotational symmetry is broken 
down to the cubic group on a finite lattice, and a mixing 
of different partial waves occurs. In moving frames, the 
symmetry is further broken down to the different sub- 
groups of the cubic group (different little groups) . Using 
symmetry arguments, it is possible to carry out a par- 
tial diagonalization of Eq. as well as to construct 
the operators that project out the spectra corresponding 
to the different irreducible representations of the little 
groups. A full-fledged analysis of the problem has been 
carried out recentl y 17311. For earlier work on the sub- 
ject, see, e.g., Refs. |29l. l34l. IHlL l56l and references 
therein. Below, we use the results of Ref. [73j |. in case of 
the meson-meson scattering in S, P, D waves, in order 
to attribute the emerging structures to the different little 
groups. 

1. The case P = (2n/L) (0,0,0) 

The partial waves S, P, D do not mix - the pertinent 
part of the matrix G is diagonal in £,£' and all entries 
with m — ml ^ mod 4 vanish. Non-vanishing ele- 
ments are: Goo.oo > Gix,n = Gio,io = Ca-^i-i , G2i,2i = 

^2-1,2-1 , G20.20 j G22,22 = G2-2.2-2 , G22,2-2 = 
G2-2.22 = (G20,20 — G21,2l)/2 • 

Below we list the irreducible representations for a given 
£ and the pertinent equations for the determination of the 
energy levels [Hj], obtained from Eq. (|33p : 



« = o, 


At: 
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— VoGoo.oo 


= 0, 


1 = 1, 


Tr ■ 


1 


— n."io,io 


= 0, 


e = 2, 


E+ : 


1 


— V^G^x^x 


= 0, 


e = 2, 


1 2 ■ 


1 


— V2G20,20 


= 0. 



Note that for more than one channel, Vi and Gt m ,t'm' m 
these equations can be considered as matrices in channel 



space. See also Sec. Ill D[ The energy levels are then given 
by the zeros of the corresponding determinants. For the 
cases of boosts with P ^ 0, one has to be slightly more 
careful as the order of the matrices matter, and the for- 
mulae of the following sections are derived for the one- 
channel case. An example of partial wave mixing with 
multiple channels is explicitly evaluated in Sec. IIV Ai 



2. The case P = (2tt/L) (0, 0, 1) 

In the unequal mass case, we restrict ourselves to S, P 
waves only. This corresponds to the partial wave mix- 
ing studied for the k(800), if* (892) system in Sec. [TTT1 
The non-zero matrix elements of G are: Goo, 00 , Gxi,xi = 
"1-1,1-1 , Gio,io j Goo, 10 = Gio.oo ■ 

The irreducible representations and the equations for 
the determination of the energy levels are given by solv- 
ing Eq. HMD: 

£ = 0,1, Ax: (l-FoGoo.ooXl-^Gxo.io) 

- VbFiG 0010 = , 
1 = 1, E: 1-^1^11,11=0. (38) 

Note that here and in the following, the Gt m .i> m > can 
be mapped to linear combinations of the Liischer zeta- 
functions or the quantities wi" m " used, e.g., in Refs. [H, 
(see Eq. ([311)). 

We shall need D-waves in the equal-mass case only. 
In particular, we will study the mixing of _D-waves with 
S'-waves for the ct(600) isoscalar in Sec. HVl For equal 
masses, Goo, 10 = and the mixing between the S- and 
P- waves vanishes. The mixing vanishes anyway for irir 
as L + I =even, and the P-wave has isospin 1. 

Additional non-zero matrix elements are: Goo. 20 = 

G20.00 j G20,20 j G22,22 = G2_2,2-2 j G21,21 = "2-1,2-1 j 
G22.2-2 = G2-2.22 ■ 

The irreducible representations and the equations for 
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the determination of the energy levels are: 

e = 0,2, A+ : (l-FoGoo,oo)(l-V 2 G 2 o, 2 o) 

- V V 2 G 0020 = , 
1 = 2, E+ : l-V 2 G 21j21 =0, 
i = 2, B+ : 1 - V 2 (G 22 , 2 2 + 622,2-2) = , 
1 = 2, B+ : 1 - V 2 {G 22 , 22 - 622,2-2) = . (39) 

3. The case P = (2ir/L) (1, 1, 0) 

For the unequal mass case, we again consider S, P 
waves only. The non-zero matrix elements of G 
are: Goo, 00 , Goo, 1-1 1 Goo, 11 = — ?Goo,i-i , Gi_i,oo = 
*Goo,i-i , G11, 00 = —Goo, 1-1 j Gio,io , Gxx,u = 
Gi_i,i_i , G11, 1-1 = — Gi-1,11 . The irreducible repre- 
sentations and the equations for the determination of the 
energy levels are: 

£ = 0,1, M ■ (l-^oGoo,oo) 

x (l-Vi((5i_ M _i-i6i M _i)) 
— VoV\ (Goo, 11 — Goo,i-i) 2 = , 

l=\, B x : 1-^1610,10 = 0, 

1=1, B 2 : l-Vi(Gi_i,i_i+tGn,i_i) = 0. (40) 

The mixing of S- and 7D-waves is again considered in the 
equal-mass case. The irreducible representations and the 
equations for the determination of the energy levels are: 

£ = 0,2, Af: (l-y G o,oo) [(1-^2620,20) 

X (1 — V2(G2-2,2-2 — G2-2,22)) 



2^2 6 2 _ 2 20 



+ 2(1 — V2620, 20)^0^2600 2-2 



Goo,2oVbV2 



G 



00,20 



X (1 — V2G2-2.2-2 + ^2G2-2,22) 









4 


- 4y 2 Goo, 2 - 2 G 2 


-2,20 — , 




1 = 


'2: 


At: 


1 


— V 2 {G 2 -X,2-1 


- i6 2 -l,2l) 


= 0, 


1 = 


2, 


Bt: 


1 


— V 2 (G 2 -X, 2 -X 


■f iG 2 _l,2l) 


= 0, 


£ = 


2. 


B+: 


1 


— V 2 (G 2 - 2 , 2 - 2 


+■ G2-2,2 2 ) = 


= 0. (41) 



4. The case P = (2tt/L) (1, 1, 1) 

For the unequal mass case, we consider S, P waves 
only. The non-zero matrix elements of G are: Goo, 00 ? 
600,1-1 = —611,00 = 2 1 / 2 (1 — «)6oo,io 7 6i_i,oo = 
— 6oo, 11 = 2 X / 2 {1 + i)6oo,io ■ 600,10 = 610,00 j 

1-1,1-1 — (j-10,10 — *jii,ii , <-ni,i-i — — (ji-i,ii j 
61-1,10 = —610,11 = 2 1 / 2 (1 — «)6ii,i-i , 6n,io = 
-610,1-1 = 2 _1 / 2 (1 + i)Gu,i_i . The irreducible rep- 
resentations and the equations for the determination of 
the energy levels are: 

£ = 0,1, Af. (l-Vi(Gi_i,i_i+2iGi_i > ii)) 

x (1 - ^0600,00) - 3^ ^iGg 0>10 = , 
1=1, E: l-1q(Gi_i,i_i-iGi_i,ii)=0. (42) 

The mixing of S- and D-waves is considered in the 
equal-mass case. The irreducible representations and the 
equations for the determination of the energy levels are: 



i = 0,2, A+ 



V V, 
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UV3G2-2.20 + 3V2G2-121 



+ (l-F 2 (G 2 _i >2 _i+2iG 2 _i j2 i)) 
x (1 - VbGoo.oo) = , 
1 = 2, E+ : (1 - F 2 6 20 ,2o) 

(1 - y 2 (G 2 _i >2 _i - iG 2 _i,2i)) 

-6F 2 2 6 2 _i,2o6 2 o,2-i =0. (43) 

As mentioned before, our results completely agree with 
those of Ref. [7311 in all partial waves. Further, note 
that in Ref. [53], the equations are derived for the cases 
P = (2n/L) (0,0,1) and P = (2^/^(1,1,0), only S and 
P waves retained. We have checked that our results ex- 
actly agree with the results of Ref. [57] for these cases 
(the representations Bx,B 2 for P = (2tt/L)(1, 1, 0) are 
denoted by B3 , B 2 in Ref. [57| ) . In Ref. [56[ , results for 
the boosts P = (2n/L)(0, 0, 1) and P = (2tt/L)(1, 1,0) 
have been obtained, and in Ref. [58| results equivalent to 
the present ones for the SD-wave case were determined. 

F. Calculation of the energy levels 

In practical calculations, we found it convenient to con- 
sider Eq. (|33|) in a rotated frame. We namely rotate the 
z direction of the coordinate system, in which (9, 0) of 
q* are measured, into the direction of the boost vector 
P. To rotate z actively into the P direction - the latter 
given by (8p, <fip) - one has different choices; we choose 
here first a rotation around the i-axis by — 9p, so that 
the new z ' vector is in the yz plane, and then a rotation 
around the original i-axis by <f>p — ir/2 so that z"||P. 
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The coordinates q * " of a vector q * , in the rotated co- 
ordinate system with z-axis z", are then given by the 
inverse of these (non-commuting) rotations. With stan- 
dard rotation matrices Ri, i = x,y, z, the new Cartesian 
coordinates read 



q* " = Rq* = [R z [(b P - -) R x (-6 P ) j 



sin (bp 



— cos 6p 







R= cos 8p cos (bp cos Op sin (bp — sin 9p 
\ sin 9p cos (bp sin#psin</)p cos Op 



.(44) 



where 

P = 0, (b P = for P = (2tt/L)(0,0,1) , 

<9 P =71-/2, p =tt/4 for P = (2tt/L) (1,1,0) , 

6»p = arctanV2, P = tt/4 for P = (2ir/L)(l, 1, 1) . 

(45) 

We would like to stress that all calculations in this pa- 
per have been carried out in the rotated frame. The 
formulae listed in the previous section merely serve to 
demonstrate the equivalence of the present approach to 
that of Ref. [z3. 

From now on we will only consider the G functions in 
the rotated coordinate system (compare to Eq. (|3Tj) ). 

47T v-^ hi f q 



£ 3 ^ PoV?° 

71 

xY; m (f"(q*))Y e , m ,(q*"(q*))I(q*) . (46) 

In particular, with the rotation of the coordinate sys- 
tem the zeros of the determinant, i.e. the levels, remain 
unchanged and the matrix of Eq. (|33[) becomes block- 
diagonal. 

For the boosts P = (2tt/X)(0, 0, 0), (2tt/X)(0, 0, 1), 
(27r/i)(0,l,0), (2tt/L)(1,0,0), (2tt/L)(0, 0, 2), 
Eqs. d23), O, (JMD remain of course valid with G R from 
Eq. (|4*6"|) , because for these boosts Gf ml , m , = Ge m ,£'m'- 
For the boost P = (2n/L)(l, 1,0), we obtain: 

£ = 0,1, A i: (l-FoGg,,oo)(l-^fo,io) 

-^i(G«, 10 ) 2 = 0, 

< = 1, l-^i(Gf 1)11 +Gf lil _ 1 ) = 0, 

£=1, B 2 : l-V 1 (G R -G R ) = 0, 



£ = 2, A+ : 1-V 2 (G 



_ s~<K 
22,2-2 ^22,22 



0. 



£ = 2, B+: l-V 2 {G R lj21 +G R >2 _ 1 ) = Q, 
£ = 2, B+ : 1 - V 2 (G R lt21 - Gf li2 _ x ) = , 
£ = 0,2, A+: l-T/ G^ 00 - a y 2 + 6F 2 2 = (47) 



where 

a = G 2 ? _2,2-2 + G 2 7 _ 2: 22 + G^^o 
+ Vb[2(G^ 02 _ 2 ) 2 + (G^ 020 ) 2 

— oo {p2-2,2-2 + G^_ 2 22 + G^q^o) 



6 = 2(G«_ 5 



K,G 



4 Gq^, 2 _ 2 Gqq i 2 q 2 ; 2 q 



{G R _ 2 .2-2 + G 2 _ 2i22 ) G 
[2G 2 | )j2 o(G ^ );2 „ 2 ) 2 



20,20 



(G 2 , _ 22 _ 2 + G 2 _ 



((G 



00,20 J 



'00,00^20,20^ 



(48) 



For the boost P = (27r/L)(l, 1,1), we obtain: 

£ = 0,1, A x ; (l-yoGS l00 )(l-ViGS,io) 
- Wi(G&, 10 ) 2 =0, 

£=1, E: l-Vi(3f M1 =0, 

^ = 0,2, A+: (l-VoG« 00 )(l-F 2 Gf 0i20 ) 
-^o^ 2 (G l )i2 o) 2 = 0, 

£ = 2, E+: (l-V 2 G R a2 ){l-V2G R 1 . 21 ) 

-V 2 2 G R j2 _ 1 G R _ ia2 = Q. (49) 

Finally, we denote that the structures of Eq. ([56]) are the 
same with rotated G R , except that the rotation renders 
many G R equal to zero, so that the block-diagonal struc- 
ture becomes immediately visible. 



G. Equations for the determination of the 
scattering phase shifts 

One of the main tasks of this study is the reconstruc- 
tion of the S-wave phase shift from lattice data taken in 
moving frames. We will assume that suitable interpolat- 
ing operators are used in lattice simulations that allow 
to associate the levels obtained to any of the symmetry 
groups discussed. A step forward in this direction has 
been given, e.g., in Refs. [HI, [73[ by constructing a basis 
of interpolating operators transforming irreducibly under 
the reduced symmetry of the moving particle in a cubic 
box. Similarly, a complete list of two-particle interpola- 
tors is quoted in Ref. [57] ] for meson momenta pi where 

Pl,P2 < 2-Ky/i/L. 

The structure of the equations for the level determi- 
nation, denoted in the previous section, lends itself to 
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a suitable strategy: the P- or D-waves can be recon- 
structed separately and serve as input to disentangle the 
S-wave from the representations A\ and A\ in which P- 
or D-waves appear mixed with the S-wave, respectively. 
For the one-channel case, a corresponding set of equa- 
tions is formulated in the following. 

Note that to disentangle partial waves one needs in 
principle eigenvalues from different irreducible represen- 
tations at exactly the same energy, which will unlikely be 
the case in an actual lattice calculation. This requires to 
make assumptions on the continuity of the amplitude so 
that one can interpolate to different energies. Strategies 
of how to do this in practice are discussed in Sec. IIVB 11 

The key point for phase extraction is that, although 
Eqs. ([SHI EH E|) and (09]) are formulated in terms of the 
divergent, and thus cut-off dependent G, the extraction 
and disentanglement of partial waves can be formulated 
entirely in terms of 



1 — G«L fi m i — Su'8 mm i ( G + 



8nE J 



GIL h m i — SpjiSm.m' Re G 



(50) 



with G R from Eq. |gBJ) and G from Eq. flU (the last line 
of Eq. (|5L)p is valid above threshold). For the quantity G 
the dependence on the cut-off cancels because it depends 
only on the difference G — G; the quantity G for the 
case 5w5 mm i = is convergent anyway: for q* much 
larger than the considered typical momenta, the sum (c.f. 
Eq. (|46p ) can be approximated by the integral and then 
Eq. p9p renders the high-momenta contributions to zero. 

Below we demonstrate in one example, how to disen- 
tangle S-waves from the P-waves. We consider the irre- 
ducible representation A\. Eqs. (EE]), (g7|) and (gSJ) all 
have the same form in this representation 

(1 - Vb<5« 00 )(1 - ViG? QtW ) - W<5« 10 ) 2 = 0. (51) 

This equation determines the A\ levels for the boosts 
P - (0, 0, 1), (1, 1, 0), (1, 1, 1). We further introduce 



Ti = 



Kt = 



Vi 



IttE 



1 - ViG p cot Si(p) - ip ' 
Vi -8nE 



i = 0,l. (52) 



l~V t ReG p cot 5i(p) 
With these definitions, Eq. (|5Tj) is rewritten as 

(l-/v Goo,oo)(l-A" 1 G 1 o,io)-A A' 1 (Goo,io) 2 = 0. (53) 

Solving this equation with respect to the S-wave Sq = S$ 
as a function of the P-wave Si = 5p, we finally obtain 

{8nE) 2 G 2 00 10 

p cot S = -8tt E Gqo.oo H — — — — ^ , (54) 



p cot Si + 8tt E Gio, 
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which can also be obtained from Eq. (|51|) upon substi- 
tution of G R by G where G is defined in Eq. (|50")l . Note 



that, although in the derivation of Eq. (|54[) we have used 
the quantities G and G, the final result depends only on 
the cut-off independent quantity G, as has to be. 

If the P-wave phase shift vanishes, Eq. simplifies 

to 



p cot<5 = — 8tt EG, 



00,00 



(55) 



i.e. the standard Liischer formula for pure S'-wave, in the 
formulation of Ref . [32J . 

Eq. (f54"} provides - in the one-channel case - the pos- 
sibility to fully correct for partial wave mixing at the en- 
ergy E of the measured level, but it requires the knowl- 
edge of the P-wave phase shift 5\ at that energy. As 
discussed above, this knowledge may come from a sep- 
arate analysis of levels with pure P-wave content, i.e. 
from other representations than A\, c.f. Eqs. ([55]) . lj4"T|) 
and (f4"9"|) . The following cases allow for this extraction 
of the P-wave: for P = (2n/L)(0, 0, 0): representation 
Tf , c.f. Eq. (E3); for P = (2tt/P)(0, 0, 1): representation 
E, c.f. Eqs. (EE}; for P = (2tt/P)(1, 1, 0): representa- 
tions Pi and B 2 , c.f. Eq. (07]); for P = (2tt/L)(1, 1, 1): 
representation E, c.f. Eq. (f4"9"|) . 

Finally, we consider the mixing of S- and P/-waves in 
case of equal-mass particles. Extraction of the D-waves 
alone proceeds analogously to the extraction of the P- 
waves in all cases except P = (2vr/L)(l, 1, 1). In this 
case, there are two independent solutions 



pcot(5 2 = —4:nE 



G*22,22 + G*21 21 



± \ 4G22, 2-1^2-1, 22 + ( ^22,22 — ^21,21 



(56) 



This is related to the fact that, in this case, there are two 
different representations E + for 1 = 1. 

Finally, the mixing of the S- and D-waves can be 
treated analogously to the mixing of the S- and P-waves, 
since the equations are always linear in cot Sq. In the case 
P = (27r/P)(l, 1,0) the final equations are a bit cumber- 
some and contain terms quadratic in cot ^2 (not explicitly 
quoted here). 

For the two-channel case, the structure of the deter- 
minant (|33[) is more complicated, because the V and G 
are matrices in channel space as discussed in Sec. Ill Dl 
There is not much point to formulate analytic formulae 
as in this section, because for a disentanglement of the 
S'-wave in a channel, one needs to know not only the P- 
or D-wave phase shift, but also inelasticities and phase 
shifts in the other channel. For S-wavc and P = 0, the re- 
construction of the amplitude in the two-channel case has 
been discussed and solved in Ref. [32J . In this study, we 
will construct synthetic data using the full two-channel 
case and reconstruct the one-channel phase shifts using 
the framework developed in this section. Below the in- 
clastic threshold, this provides a very good approxima- 
tion whose validity will be discussed. 



10 




700 800 



900 1000 
E [MeV] 



1100 



FIG. 1: Solution for the isospin I = 1/2, strangeness S = —1 
meson- meson interaction (solid lines) of Ref. [iff]. Partial 
wave data for 5®/ 2 : circles (75|, triangles up: average as de- 
fined in ref. [Zf|. For 5\, 2 : triangles up fz?}, circles 



III. THE k(800)/K*(892) SYSTEM 
A. Level spectrum 
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FIG. 2: Spectrum of the S- and P-wave system in coupled- 
channel scattering {irK, r/K) for P — (0,0,0). Solid (dashed) 
lines: S'-wave (P-wave). The fine dotted lines show the non- 
interacting levels. 



serve an S-wave level close to the wK threshold as well 
as avoided level crossing at the r\K threshold for the S- 
wave. These features have been discussed extensively in 
Refs. [32l. [45l] . In particular, the avoided level crossing is 
a signal of the r]K threshold (c.f. fine dashed horizontal 
line) and not of the k(800). 



To study the mixing of partial waves and apply the 
framework developed in previous sections, the case of 
coupled-channel irK, r\K scattering in a moving frame 
is considered. In isospin I = 1/2, strangeness S = — 1, 
we take account of the S-wave with the k(800) resonance, 
the P-wave with the K*(892), and neglect all higher par- 
tial waves. S-wave irK scattering and the n resonance 
have been addressed in recent lattice calculations, e.g. in 
Refs. Hi,|27|]. 

For the hadronic amplitude in the infinite volume, we 
use the solution obtained in Ref. [45| by fitting the low 
energy constants L\ to L% to strangeness S = 0,-1 par- 
tial wave data in S- and P-waves. That solution is ob- 
tained using the inverse amplitude method of Ref. [62| . 
In Eq. (|34[) the connection to the present framework is 
quoted. The phase shifts are shown in Fig. [TJ For the 
S'-wave we observe a pronounced effect due to the r\K 
threshold which shows that this channel can have notice- 
able influence. For the P-wave, the r/K channel plays 
almost no role. 

Using the equations of Sec. lIIEI or lllFl the lattice spec- 
trum for the different boosts can be predicted as a func- 
tion of the box length L. For P = 0, the result is shown 
in Fig. [5J Here, S- and P-waves do not mix. We ob- 



As we choose a finite boost P 7^ fj, the level spec- 
trum becomes more complex. For P = (2ir/L)(0, 0, 1), 
P = (2tt/L)(1,1,0), and P = (2tt/P)(1, 1, 1) the result- 
ing levels as a function of L are shown in Figs. [31 [3J and 
respectively. 

As the boost increases, the (non-interacting) levels are 
in general moved upwards in energy. Close to the low- 
est boosted non-interacting level, i.e. the boosted 7rA" 
threshold, one always finds a level in which S- and P- 
waves are mixed (solid lines, representation Ai) . At 
higher energies, apart from levels of this kind, also lev- 
els from other representations start to appear containing 
only P-wave without S'-wave admixture and allowing for 
a separate analysis of the P-wave. 



The levels of the representation A\ are determined from Eq. J5l}. 
The first term of Eq. II51I I suggests the occurrence of two levels 
that undergo a mixing through the second term ~ V0V1. How- 
ever, as e.g. Fig. [3]shows, there is only one level close to the first 
non-interacting energy, and only one close to the second non- 
interacting energy. Indeed, it can be shown that this must be 
the case, by using the known angular structure of the G( m y m i . 
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FIG. 3: Spectrum of the S- and P-wave system in coupled- 
channel scattering (ttK, i]K) for P = (2n/L)(0, 0, 1). Solid 
lines: levels with SP-wave mixing from the irreducible rep- 
resentation A\. Dash-dotted (red) lines: levels from only P- 
wave, irreducible representation E. The fine dotted lines show 
the boosted non-interacting levels. 



1100 



1000 - 



900 



800 - 



700 




FIG. 4: Spectrum of the S- and P-wave system in coupled- 
channel scattering (ttK, i]K) for P = (2n/L)(l, 1, 0). Solid 
lines: levels with S P-wave mixing from the irreducible rep- 
resentation Ai. Dash-dotted (red) lines: levels from only P- 
wave, irreducible representations B\ and E>2. The fine dotted 
lines show the boosted non-interacting levels. 



B. Leading behavior of the level shifts 

The levels from the A\ representation get shifted once 
the mixing is taken into account. To understand this 
effect, we consider mixing in the limit of small phase 
shifts. In the absence of a P-wave (V\ = 0), Eq. (|?T|) 




UM„ j 



FIG. 5: Spectrum of the S- and P-wave system in coupled- 
channel scattering (ttK, rjK) for P = (27r/L)(l, 1, 1). Solid 
lines: levels with SP-wave mixing from the irreducible rep- 
resentation Ai. Dash-dotted (red) lines: levels from only P- 
wave, irreducible representation E. The fine dotted lines show 
the boosted non-interacting levels. 



reduces to 

1-V S G^ 00 (E = E S ) = (57) 
which determines the position E$- The level shift is 



AE = E S p- E s 



(58) 



with Esp the position of the level with mixing, given by 
the solution of Eq. (|5ip . The shift can be approximately 
calculated by Taylor expanding Eq. ([ST]) around E = Es- 
It is straightforward to show that then approximately 



AE ~ S, 



8irE s (G^o.io) 2 



P dG*/dE 



(59) 



where 5p is the P-wave phase shift, p the three- 
momentum of the 7r and K in the two-particle rest frame, 
and G defined in Eq. (|5U|) . The actual shift AE is shown 
with the solid lines in Fig. [HI the result from Eq. ([59")) is 
indicated with the dashed lines. Indeed, for the lowest 
level where Si is small Eq. (|59p provides a good approx- 
imation. 

To understand the structure of the shift, we can expand 
G around the pole at E = Eq where Eq is the boosted 
non-interacting energy, 



a_i h 



(n=3) 



E-E 



+ Rtr, 



(60) 
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FIG. 6: Level shift through partial wave mixing for P — 
(2ir/L) (0,0,1). Below: for the lowest level from Fig. [3] 
Above: for the following level. The solid lines show the ac- 
tual level shift, the dashed lines the result of Eq. (|59p and the 
dotted lines the leading behavior from Eq. I|62|). 



where h 



(<>) 



(n=0) 



contains the angular structure, 



a_i h\ m e / m , is the residue, and Re m ,£'m' is regular 
in the vicinity of Eq. The residue is readily evaluated, 



1 



a_i 



1 



L 3 2W1W2 E + OJi + Ll>2 



(61) 



where the energies u>i are evaluated with the boosted 
vectors q* = q*(q = ^(0,0,0)) for the lowest level 
and q*(q = ^(0,0,1)) for the following one. Using 
uii + L02 — E and substituting Eq. ([501) in Eq. ([59]) we 
obtain 



AE 



6irEg Sp 
L 3 pu>iui2 



(62) 



for the level shift. The result of this approximation is 
shown with the dotted lines in Fig. [SJ This equation 
shows that the shift is to leading order in E proportional 
to 5p and to 1/L~ 3 . This L~ 3 behavior is similar to the 
one of the scattering length [H, HH . While Eq. (JB2J) is 
useful to understand the qualitative behavior of the level 
shift, its quantitative use is limited. 



C. Disentangling the rc(800)/K* (892) system 

The S-wave can be extracted from lattice data of 
the representations (boost P = 0) and A\ (boosts 
P = (27r/i)(0,0,l), (2w/L) (0,0,1), (27r/i)(l,l,0), 



(2tt/L)(1,1,1)). While for P = the S-wave does not 
mix with the P-wave, for the higher boosts mixing oc- 
curs in the A\ representation, c.f. Eqs. ([571 1551 Wf \ |4"5 |) . 
We have derived Eq. ([54"|) that allows to disentangle the 
S-wave for the one-channel problem provided that the 
P-wave is known from an analysis of other levels, as dis- 
cussed in Sec. Ill Gl 

In this section, we apply Eq. (|54|) to the levels cor- 
responding to A\ shown in Figs. [3J to [51 It should be 
stressed that those levels have been obtained with the full 
two-channel formalism as described in Sec. Ill Dl based on 
the global fit of low energy constants to S- and P-wave 
partial wave data as described in Sec. IIII M In contrast, 
we will extract the phase shift from these levels using the 
one-channel equation ([54")) . Below the inelastic thresh- 
old, given by the r\K channel, this is expected to be a 
good approximation. Indeed, in Refs. [3l|, HH, the so- 
called pseudo-phase, i.e. the phase extracted with a one- 
channel formalism from a two-channel problem, provides 
an excellent approximation to the actual phase up to en- 
ergies close below the inelastic threshold. 

In Fig. [7J the actual irK S'-wave phase shift, from the 
full coupled-channel system, is shown with the solid (or- 
ange) lines, identical to the corresponding curve in Fig.Q] 
Consider the first two A\ levels, tied to the first three 
boosts [solid (black) lines in Figs. [3J [3J and[5J. In those 
levels, the S'-wave mixes with the P-wave, and Eq. (pylj) 
provides the possibility to disentangle the S-wave by us- 
ing the known P-wave shown in Fig. [TJ The result is 
shown with the long-dashed lines in Fig.[7J In the figure, 
we also indicate to which box size L the extracted phase 
shift corresponds (c.f. again Figs.[3J[IJ and[5jto see the 
connection of energies and box size, given by the levels). 
As a test of the formalism, we have shown that once lev- 
els are generated from the hadronic model in a reduction 
to one channel, Eq. (|54|) ensures the exact reconstruction 
of the phase, as must be. 

In general, the agreement with the original phase shift 
is excellent. Only at higher energies there are small de- 
viations, coming from the more and more important rjK 
channel. It should be noted that the effect on the recon- 
struction of the phase is small but the reconstruction of 
the pole position might be affected. This has been shown 
in Ref. ^ for the k(800) and P = 0. 

Second, we quantify the effect if partial wave mixing is 
neglected (in the one-channel extraction scheme we are 
are using). In that case, Eq. (|54|) reduces to Eq. (|55|) 
which is the ordinary Luschcr equation for a boosted 
S-wave system. The results are shown with the dash- 
dotted lines in Fig.[7J For very large L, we indeed observe 
that the extracted and the actual phase are similar. In 
other words, the P-wave decouples from the S-wave in 
the infinite volume limit as must be. However, even for 
L > 3M" 1 we already observe large deviations and for 
smaller box sizes than L <~ 2M" 1 , a reliable extraction 
of the phase, let alone the pole position of the k(800), is 
not possible any more. 

In summary, it is crucial to take effects from partial 
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FIG. 7: Actual phase shifts [solid (orange) lines], identical 
to the result shown in Fig. [T] and the extraction from the 
first and second level, using Eq. (|54[1 for the disentanglement 
(dashed lines) or Eq. (|55p in which the partial wave mixing is 
neglected (dash-dotted lines). Results are shown for the first 
three boosts ~ (0,0, 1), (1, 1,0), and (1, 1, 1). In the figure, 
also the corresponding values of the box size L are indicated. 
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FIG. 8: Levels of the (7 = L = S = 0) meson-meson system 
[ct(600) and /o(980)] for different boosts. The symbols at 
L = 2.5Af~ 1 and L = 3.5M" 1 indicate the 20 synthetic data 
points taken for the reconstruction of phase and <r(600) pole 
position [results in Figs. 1101 and 111]. 



wave mixing into account while effects from inelastic 
channels can be neglected to a good accuracy as long as 
one considers energies only below the inelastic threshold. 
As a side remark, we would like to stress that in different 
physical contexts it might be crucial to use a two-channel 
extraction scheme, in particular if thresholds are close to 
each other as in case of the A(1405) [13, EH El] , or if one 
wants to extract phases and poles close to thresholds as 
in case of the / (980) [Hj]. 



IV. THE ct(600) IN A MOVING FRAME 
A. Level spectrum and resonance extraction 

The first few levels for the I = L = S = quan- 
tum numbers are shown in Fig. [5] for the first five boosts. 
To obtain these levels, the S'-wave two-channel potential 
(irir, KK) from the inverse amplitude method, Vg = 
V 1AM , given in Eq. (JMI, is used in Eq. (33]) to deter- 
mine the levels. For y lAM , the corresponding solution of 
Ref. [45| for the underlying hadronic interaction is used. 

To quantify the expected error from partial wave mix- 
ing, we have calculated the level shift due to the S'-wave 
mixing with the small isospin zero D-wave. Note that for 
the considered equal- mass case (7T7t and KK), there is no 
SP- wave mixing even for P ^ 0. Anyway, the 7T7t sys- 



14 



tern in P-wave has isospin one. For the SD-w&ve mixing, 
the phenomenological parameterization of the P-wave of 
Ref. [l^] is used that serves to construct a P-wavc to D- 
wave transition V D 2 ' 2 ^ as described in Appendix[A] We do 
not include a KK channel in P-wavc as there is no phe- 
nomenological reason for it and we are far below the KK 
threshold, anyway. We have, for the representation, 
the levels given by Eq. (|33p where 
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for the boosts P = (2tt/£)(0, 0, 1), (2tt/£)(1, 1, 1), and 
(2tt/L)(0, 0, 2). For the boost P = (2tt/L)(1, 1, 0), 
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where the index (1) labels the A'-FT channel and (2) the 
7T7T channel. The resulting SP-wave mixed levels of the 
A^ representation are shown with the (red) solid lines in 
Fig. [5J together with the S'-wave levels without mixing 
from Fig. [SJ As the figure shows, the level shift from 
mixing, AE can reach up to 10 MeV at L = 2M~ X , but 
stays below 3 MeV for L = 2.5A/„T 1 and larger. The shift 
AE increases not only for smaller L, but also for higher 
energies where 8d is larger. Eq. (|62p indeed shows that 
this is the expected behavior. 

Note that the level shift from the .SP-mixing is much 
smaller than the one of the irK, r]K system discussed in 
Sec. because the irir D-wave is much smaller than the 
irK P-wave in which the if* (892) resides. 



B. Extraction of the <x(600) 

As the level shift from mixing is smaller than 3 MeV 
for L = 2.5M" 1 and larger, we concentrate on these box 
sizes and can neglect the mixing; for the reconstruction of 




MM, 1 



FIG. 9: Levels from the representation as in Fig. [8] In 
addition, the (red) solid lines show the levels mixed with D- 



the cr(600) we will assume 10 MeV errors on the synthetic 
data so that this assumption is safe. In this section, we 
concentrate on the extraction of the cr(600) resonance 
and do not analyze the /o(980). The latter resonance 
would require a two-channel extraction scheme that has 
already been discussed in Ref. [32j , and the extension to 
moving frames is in principle straightforward. 

As Fig.[5]shows, with larger boosts the first and second 
level from the A^ representation move towards higher en- 
ergies. In this way one can cover the entire energy region 
from threshold up to E = 700 MeV with the boosted 
first two levels. This demonstrates the advantage lying 
in the use of moving frames: in the conventional Luscher 
approach, L is varied at P — 0. This means that, first, 
many different lattice setups have to be calculated. Sec- 
ond, for the box sizes considered 2 < L < 3.6M" 1 , the 
energy region around E ~ 400 — 500 MeV, i.e. precisely 
where the real part of the pole position of the o (600) is 
located [8(|, is not covered by any level, making the ex- 
traction of the <r(600) generically more difficult, as has 
also been noted in Ref. [451 ] . 

In the absence of actual lattice data, we generate 20 
synthetic data points from the levels of Fig. [5] at two 
values of L, L = 2.5M~ 1 and L = 3.5M" 1 , as indicated 
with the dots in the figure. A 10 MeV error is assigned 
to each data point. With such-defined data, the task is 
to reconstruct the 7T7t phase shifts and the c(600) pole. 

For this, we use the parameterization of the one- 
channel potential from Ref. [45j , 



V 



in 



at 



Vf* =a + b(s- s ) + c(s - soY + d{s - s ) d (65) 

with V2 = Vlo the fixed LO term of the chiral expansion. 
In other words, we take the form of the inverse amplitude 
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(c.f. Eq. ([Ml) ), leaving the LO term V2 as given by chiral 
symmetry, and expand V4 in powers of s. As expansion 
point we choose s = (400 MeV) 2 . 

The choice of this potential and its advantages have 
been extensively discussed in Ref. [45[. We denote here 
that the explicit inclusion of the model- independent, 
well-known lowest order term V2 greatly helps stabilizing 
the extraction. The higher powers of s account for cor- 
rections from the next-to- leading and higher order terms. 
Note that the polynomial NLO contributions [12] can be 
approximately taken account of, because the Mandelstam 
variable t and u can be expanded in s. The effect from 
the KK channel and its branch point at E = 2Mk, as 
well as the left-hand cut, is also absorbed in this expan- 
sion. 

One could in principle include these non-analyticities 
explicitly in the fit potential. However, they are not well 
fixed because they lie much higher or much lower in en- 
ergy. Given lattice data in a relatively narrow window in 
energy, no improvement is expected but instead large cor- 
relations of the corresponding new parameters will arise. 
The expansion of these effects in a power series in s, as 
provided in Eq. (|65|) . allows for a systematic improve- 
ment and provides a set of parameters with relatively 
small correlations. 

One should bear in mind, though, that this is an ap- 
proximative procedure. Other than in Ref. [32j|, where 
the analytic form of the fit potential comprised the as- 
sumed hadronic interaction - this was possible because 
the model interaction was from lowest order only - the 
potential from Eq. (|rj5]) can only approximatively absorb 
the discussed effects. The strategy is then to perform dif- 
ferent fits with increasing powers of s until convergence 
is observed as discussed in Ref. [45[ in detail. 



1. Extraction strategies with partial wave mixing 

For the tttt 5-wavc sufficiently below the KK thresh- 
old, we can neglect the partial wave mixing as discussed 
before, but we give an outlook how to proceed in case 
it cannot be neglected as for the k/K* system discussed 
in Sec. IIIII On one hand, a separate extraction of P- or 
D-waves is possible from the structure of the irreducible 
representations B\, P2, E (P-wave) and A% , B^ , B%, 
E + (_D-wave) as discussed following Eq. ([55| . see also 
Ref. [57] (we assume there is no mixing between P- and 
D-wave). However, those phase shifts can only be ex- 
tracted at scattering energies different from where the 
levels are situated, from which the 5*-wave can be re- 
constructed via Eq. ([51j) in the one-channel case. For ref- 
erence, see Figs.[3]to[5J It is then necessary to make some 
minimal assumptions on the P- or P-wave phase shifts, 
such that the underlying potential can be expanded in 
energy as discussed in Sec. HVBl c.f. Eq. (frT5) . This has 
also been recognized in Ref. [57] . 

Instead of determining first the higher partial wave and 
then the S-wave, it might be advantageous to simultane- 



ously fit the different levels with Vs and Vp (or Vs and 
Vd), both of them expanded in energy as in Eq. frlS]) . 
Such a procedure would be in analogy with the two- 
channel extraction scheme developed in Ref. and is 
expected to be more efficient: the simultaneous fit of pure 
P-wave and S P-wave levels leads to smaller uncertainties 
in the infinite-volume limit than the two-step procedure 
discussed before. 

A generalization to the multiple partial wave, multiple 
channel situation, using the corresponding matrices in 
channel- and partial wave-space as those of Eqs. ([551 l3"o| . 
is in principle straightforward but will require very high 
precision lattice data; see Ref. [32| where this issue is 
discussed for the two-channel case. 

As performed in Sec. IIV CI for the er(600), there is 
no particular problem to generate pseudo-data for the 
k/ K*-system and disentangle phase shifts following the 
strategy formulated here. However, to not overload this 
study, we concentrate on the somewhat simpler case of 
the cr(600) in the next section. 



C. Results 

The analysis of the <r(600) proceeds as described fol- 
lowing Eq. (|65[) and in Ref. [45f . with the fit potential 
V Rt from Eq. IJB5J). The fits are labeled according to 
the powers of s used in V . Once a minimum is found, 
parameter errors are determined. The parameter error 
for a parameter a is defined by the range of a in which 
X 2 < Xbcst + 1) under the constraint that all other pa- 
rameters are optimized. Then, within the errors, random 
parameter sets are generated and only those sets kept for 
which x 2 < Xbcst + 1- For each of these sets, phase shift 
and pole position are calculated. The resulting bands 
and areas, for phase and pole position, respectively, are 
shown in Figs. [Till and [TT1 Note that the uncertainty area 
for the pole position in the (s°) fit, shown in Fig. [Til 
shrinks to a line. For the best x 2 1 P°^ e positions are in- 
dicated with symbols in Fig. [TTJ The actual phase shift 
and pole position, derived from the hadronic interaction 
that was used to generate the synthetic data, are also 
indicated in the figures. 

As visible for phase shifts and pole positions, fits 
with larger number of free parameters result (trivially) 
in larger uncertainties, and, of course, in a better \ ■ 
For the central values of the pole positions, we observe 
that with an increasing number of parameters, the ac- 
tual and the fitted pole positions get closer, but even the 
4-parameter fit does not perfectly match the actual pole 
position although it has the best x 2 of all fits. As dis- 
cussed in Sec. 1111 there are terms of higher order in s N , 
N > 3, in the original potential that cause this small but 
finite discrepancy. In the fit of actual lattice data, one 
will have the same effect, of course. 

As discussed in Ref. [45|| the effect of the heavier chan- 
nel - given by KK in the present case - usually can be 
well absorbed in the coefficients of the expansion of the 
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FIG. 10: Extracted phase shifts (bands) using synthetic data points as indicated in Fig. [8] with a 10 MeV error for each data 
point. The actual phase shift is shown with the dashed lines. 
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FIG. 11: Pole of the <r(600) extracted from the synthetic data 
shown in Fig. [S] The actual pole position is indicated with the 
large square. Extracted pole positions for the different fits: 
triangle down: (s°) fit; triangle left: (s^s 1 ) fit; diamond: 
(s , s 1 , s 2 ) fit; circle: (s°, a 1 , s 2 , s 3 ) fit. Also, uncertainties are 
shown (shaded areas). The dashed line shows the uncertainty 
from Ref. [iEl ] , where - instead of two volumes and five boosts 
as done here - 6 different volumes had to be considered. 



fit potential, but as the cr(600) is very broad, small re- 
maining discrepancies become large far in the complex 
plane. The same behavior was found in Ref. [4oT |. where 
L was varied to extract the er(600) pole. The situa- 
tion cannot be improved by explicitly including the KK 
channel in the extraction process, in contrast to the case 
of the /o(980) where this is possible and necessary |32j]. 
Here, i.e. still far below the KK threshold, the channel 



transitions VL 



and V RK _ 



■kk are very weakly con- 



strained. The large number of new free parameters, tied 
to these additional transitions, would immediately lead 
to drastically increased uncertainties on the observables 
and large parameter correlations as has been tested. 

We observe, in any case, that using synthetic data 
from values of L smaller than 2.5M~ , i.e. higher en- 
ergies E, immediately helps to narrow down the uncer- 
tainties of the phase shifts at higher energies, shown in 
Fig. [TO] However, in that case the central value of the 
pole position in the (s°, s 1 , s 2 , s 3 ) fit, shown with the cir- 
cle in Fig. [TT] starts to deviate considerably from the 
actual pole position (empty square), which is a sign of 
the increasing effect of the KK channel on the lattice 
data closest to the KK threshold. Choosing box sizes 
of L = 2.5AI~ 1 and L = 3.5M" 1 as done here provides, 
thus, a good compromise. These values of L, for which 
partial wave mixing can be safely neglected, is promising 
for the setup of an actual lattice simulation to extract 
the a (600). 

Finally, we would like to compare the present results to 
those of Ref. [45| where - instead of two required volumes 
as is the case here - 6 different volumes had to be con- 
sidered, implying a much increased numerical effort for 
actual lattice calculations. The uncertainty of the <j(600) 
pole position, coming from synthetic data with the same 
10 MeV error as used here, is shown with the dashed line 
in Fig. [TTJ The extension is of similar size as the one of 
the (s° , s 1 , s 2 , s 3 ) fit performed here. 

Thus, with only two different volumes and using data 
from five boosts for each, one can expect results of simi- 
lar precision than from six different values of L without 
boost. This demonstrates that the proposed extraction 
method can be quite effective in the analysis of actual 
lattice data. 
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FIG. 12: Uncertainties of the <r(600) pole position from the 
statistical uncertainty of the data. As indicated with the ar- 
rows, the uncertainties increase compared to those of Fig. 1111 
The cases of the (s°, s , s 2 ) and the (s°, s , s 2 , s 3 ) fits are 
shown. 



D. Statistical error 

One should note that for the generation of the syn- 
thetic data we have so far only assigned an error, but not 
allowed the statistical fluctuation of the centroids of the 
error bars. In other words, we have assumed that the en- 
tire error is systematic. Here, we allow, in addition to the 
10 MeV error, a statistical fluctuation of the centroids by 
5 MeV, as has been done in Ref. [12] for the extraction 
of the Jo (980). The values of 10 MeV for the error bar 
and 5 MeV for its fluctuation are chosen to demonstrate 
the effect; when it comes to the analysis of actual lat- 
tice data, these values and in particular the fraction of 
the statistical in the total error have to be adapted, of 
course. Using the same procedure as in Ref. [3^ | , we can 
estimate the resulting uncertainties, shown in Fig. [T^] for 
two fits. As the figure shows, the uncertainties increase 
by around one third for the chosen values. 



V. SUMMARY 

The present study provides a formulation for the scat- 
tering of two particles confined in a finite box with total 
nonzero momentum, adapted to the chiral unitary frame- 
work. The idea is based on extending previously known 
techniques for zero momentum, discretizing the energy 
levels by imposing the boundary conditions in the mov- 
ing frame. 

Given a hadronic interaction, levels for the first five 
boosts P = (27r/£) (0,0,0) to (0,0,2) can be predicted 
and attributed to the subgroups of cubic symmetry Em- 



ploying coupled-channel unitarized chiral perturbation 
theory including NLO terms, we derive the levels for the 
mixed-partial wave system with 1=1/2, S = —1 and 
L = 0, 1 [k(800) and IT* (892), respectively] as well as for 
the scalar sector with 7 = 0,5 = and L = 0, 2 where 
the ct(600) resides. 

We demonstrate for the k(800)/A'*(892) system that 
partial wave mixing is a very large effect for realistic box 
sizes and needs to be taken into account. To disentan- 
gle the iS-wave from P- or D-wave, we derive a set of 
equations in the one-channel formalism that are shown 
to be very precise as long as one stays below the inelastic 
thresholds. 

Furthermore, we present a scheme in which the 
hadronic interaction is expanded in energy to allow for 
the extraction of the infinite volume limit, simultaneously 
fitting levels for different boosts and at different energies. 
The model-independent information from the lowest or- 
der in the chiral expansion is kept explicitly in this ex- 
pansion, greatly stabilizing the fit to lattice data. Such 
statistical analyses can be used for actual lattice data, 
or, as done here, serve to determine promising lattice se- 
tups and the accuracy of lattice data to allow for reliable 
resonance extraction. 

The method is tested for the example of the <r(600). 
First, we show that for L > 2.5 M~ x effects from SD- 
wave mixing can be neglected. Second, synthetic lattice- 
data are produced and analyzed. We find that with only 
two different box sizes one can expect a similar precision 
on the <t(600) pole position as by varying as much as 6 
different box sizes at zero total momentum. 

Using information from moving frames is, thus, indeed 
rewarding since, with only a few different box sizes, phase 
shifts and resonance parameters of excited mesons can be 
determined. 
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Appendix A: Parameterization of the D-wave from 
Ref. (79[ 

For convenience we quote the parameterization of the 
isospin zero D-wave from Ref. [79(. To construct the 
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(2 2) 

tttt D-wavc to D-wavc transition Vp , we express it in 
terms of cot given in Eq. (A9) of Ref. (79j . 



RcG 



p cot 6 



D 



8irE 



(Al) 



where M h = 1275A MeV, B = 12.40, B t = 10.06, and 
p is the relative momentum of the pions. 



valid up to the KK threshold (E < 2M K ). Here, 

E 
V 



jt S D = — {Ml - s) Ml [B + B 1 w{s)] 



iu(s) = V ,^o = 1050 MeV (A2) 
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